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Abstract: 

The rate of quantum ergodicity is studied for three strongly chaotic (Anosov) systems. The 
quantal eigenfunctions on a compact Riemannian surface of genus g = 2 and of two triangular 
billiards on a surface of constant negative curvature are investigated. One of the triangular 
billiards belongs to the class of arithmetic systems. There are no peculiarities observed in the 
arithmetic system concerning the rate of quantum ergodicity. This contrasts to the peculiar 
behaviour with respect to the statistical properties of the quantal levels. It is demonstrated 
that the rate of quantum ergodicity in the three considered systems fits well with the known 
upper and lower bounds. Furthermore, Sarnak's conjecture about quantum unique ergodicity 
for hyperbolic surfaces is confirmed numerically in these three systems. 
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I Introduction 



In classical mechanics the behaviour of systems can be classified in the range from integrable 
towards strongly chaotic. Thereby one of the lowest suppositions of chaos is given by ergodicity 
being followed by the mixing property and finally by the Anosov property. Thus as a first 
step towards chaos it suggests itself to study the quantum mechanical analogue of the classical 
ergodicity theorem (see e.g. ]![]). 

Here, we restrict us to bounded Hamiltonian systems having a purely discrete quantal eigen- 
value spectrum {E n }, n £ N, with corresponding eigenstates \P n which are assumed to be 
orthonormalized with respect to the L2 norm. Especially, we are concerned with the eigen- 
states of the Laplace-Beltrami operator A on compact Riemannian manifolds M with constant 
negative curvature, i. e. 

- Atf„ = E n m n . (1) 

In quantum mechanics the following measure was introduced for compact Riemannian manifolds 
M 1,1 

S k (E;A) := £ | ( A* n , tt n ) - a A \ k , k > , (2) 

where N(E) := #{E n \E n < E} is the staircase function counting the number of the quantal 
levels E n below the energy E. The operator A is a pseudo-differential operator of 0-th order, 
i. e. roughly speaking, it is an operator containing no derivatives. For an introduction to the 
theory of pseudo-differential operators, see j|, ||. The average a a is defined as 



a A ■■= TTBl / cr A dT , (3) 
<S*M 



vol(S*M) 



where S*A4 is the unit cotangent bundle, i. e. the energy shell in phase space. The Liouville 
measure on S*M is denoted by dr, and a a is the principal symbol of A, which can be considered 
as the classical analogue of A. If the flow on S*M. is ergodic it is proven || ^, ||, |9| that 

S k (E;A) -> for E^ 00 . (4) 

Thus the vanishing of (|2|) in the limit E — > 00 can be considered as a quantum analogue of 
ergodicity. 

The definition of quantum ergodicity as (Q) is equivalent to another formulation. For bounded 
sequences {a n }, a n G C, one has the following lemma The condition 

lim — I a n \ p = for any p > 

n<N 

is equivalent to 

lim a n = , 
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where {rij} is a subsequence of N of density 1. A subsequence rij, n\ < nz < . . ., is defined to 
be of density d (0 < d < 1) if 

ihn iphi^a = d . 

jv^oo #{re|n < N} 

Thus the condition (||) of quantum ergodicity reads as 

(A^ nj ,^ nj ) -> cta for £ nj -» oo (5) 

for a subsequence {n?} of density one. The kind of the decay of @ for E — > oo and @ is coined 
as the rate of quantum ergodicity. Quantum unique ergodicity happens to be if the condition 
dH) holds for the complete sequence {re} 

(A^ n ,^ n ) -> cr A for £ n oo . (6) 

Attention has to be paid for interpreting the condition (|j) as a measure for quantum ergod- 
icity, since one cannot conclude from (Q) that the classical system is ergodic. An instructive 
example is the sphere S 2 , which can be considered as the counterpart of the surfaces studied 
in this paper since the sphere S 2 has constant positive curvature. The usual basis of spherical 
harmonics {Y™} is indeed not quantum ergodic as becomes clear by considering subsequences of 
a fixed rational m/l. However, for almost all orthogonal bases the eigenfunctions of the sphere 
S 2 are quantum ergodic [11]. An example of a subsequence of density zero violating (||) is pro- 



vided by the subsequence {Y[ 1 } which is very strongly localized at the equator. If one considers, 
however, the subsequence {5^ m } with m/l > c, < c < 1, one obtains a subsequence of density 
> which is concentrated in the neighbourhood of the equator given by sin# > c JP^ ]. 
In the case of compact Riemannian manifolds of negative curvature the best known upper 
bound of Sk(E; A) is given by g 

S k (E;A) = o(QogE)-t) . (7) 



By considering the sum over eigenstates in eq.@ one avoids possible sparse subsequences {rij} 
of zero density, for which 

| (A* nj ,* nj .) -a A \ = , 
which might correspond to eigenstates scarred by short classical periodic orbits. However, this 



does not occur for the hyperbolic octagons studied in jl3[. Here is the Hardy-Littlewood- 
lower-bound symbol being defined as 

f(x) = Q(g(x)) for x — > oo if f(x) / o(g(x)) . 

If the Liouville average of the sub-principal symbol <7 su b(^4 ) does not vanish, then a lower bound 
is given by || 

| (A* nj *„) - a A | = n[E~^) . (8) 
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Addressing the question of quantum ergodicity in individual eigenstates, a conjecture is put 



forward by Sarnak [14] which states that for hyperbolic surfaces in which the curvature is 
negative 

(A¥ n ,tf n ) -o A = O (e;} +£ \ for all e>0 . (9) 



The main topic in this paper deals with the rate of the decay of S\{E;A). An extended 
analysis including higher Sk(E;A), k > 1, and other statistics of the wave functions of the 
hyperbolic billiards considered in this paper can be found in [fL5f | . A detailed analysis concerning 
the behaviour of Si(E;A) of Euclidean billiards, i.e. the stadium billiard, the cardioid billiard 
and the so-called cosine billiard, can be found in . 

II The chaotic systems 

Surfaces of constant negative Gaussian curvature constitute a simple model for chaos since the 
classical dynamics on these surfaces possess all of the characteristic properties of instability: 
ergodicity, mixing and the Anosov property. As a model for this surface we use the Poincare 
disc T> which consists of the interior of the unit circle in the complex z-plane (z = x\ + ix-z) 
endowed with the hyperbolic metric 

4 

9l 3 = {l- x \-xl)2 5i i ' i '-? = 1 ' 2 ( 10 ) 

corresponding to constant negative Gaussian curvature K = — 1. This fixes the length scale. 
The systems which are studied in this paper are defined on this hyperbolic surface. 

The classical motion (geodesic flow) is determined by the Hamiltonian H = nzzPig^Pj, 
Pi = m gij dx 3 / dt. The geodesies are circles intersecting the boundary of the Poincare disc T> 
perpendicularly. 

The quantum mechanical system is governed by the Schrodinger equation 



A* n (*) = E n * n (z) , A = i(l-x?-*2) 2 (J^ + JL) t (ii) 



where A denotes the non-Euclidean Laplacian corresponding to the hyperbolic metric (1C). 
Here and in the following units h = 2m = 1 are used. The wave functions ^f n (z) have to obey 
appropriate boundary conditions, see below. 

Since we are interested in possible differences in the rate of quantum ergodicity between 
non-arithmetic and arithmetic systems, we choose the asymmetric octagon shown in figure |l] as 
the non-arithmetic system, and the so-called regular octagon for the arithmetic system. Both 
are Riemannian surfaces of genus two. (For more details about these Riemannian surfaces, 
see e.g. 17, [l8|].) However, both systems decompose into irreducible symmetry representations 
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such that we consider only one symmetry class. In the case of the non- arithmetic octagon which 
only possesses one involution symmetry the positive symmetry class is chosen, i.e. *&(—z) = 
+^(z). The regular octagon which has the highest symmetry possible for a surface of genus 
2, decomposes into 96 symmetry classes from which the one is selected which corresponds to 
a hyperbolic triangle with angles a = 7r/3, = ir/8 and 7 = tt/2. The sides a, b and c being 
opposite to their corresponding angles are endowed with the boundary conditions Dirichlet, 
Neumann and Neumann, respectively. Since these two desymmetrized systems differ in some 
aspects, we incorporate in our study a further triangular billiard which differs only in the imposed 
boundary conditions from the former. For the sides a, b and c the boundary conditions Neumann, 
Neumann and Dirichlet are chosen, respectively, which do not provide a subspectrum of the 
regular octagon. The main point is that this combination of boundary conditions leads to a 
non-arithmetic system. For a detailed investigation of these two hyperbolic triangular billiards 
with respect to the quantal level statistics, see [19|. 

Thus our study is based on one parity class of a non-arithmetic Riemannian surface and on 
one hyperbolic triangle endowed with two different boundary conditions from which only one is 
arithmetic. In the following the non-arithmetic billiard is called billiard A and the other billiard 
B, see figure [j]. Using the boundary element method [ 20 1 the first 3000 wave functions of the 
positive symmetry class of the asymmetric octagon are computed. For billiard A and B the first 
2092 and 2099 wave functions are determined, respectively. The following statistical analysis is 
based on this large set of wave functions. 



Ill On the behaviour of Si(E;A) in the configuration space 

The pseudo-differential operator of 0-th order A is chosen to be the multiplication operator with 
the characteristic function of a subset X C A4 on the Poincare disc T>, i.e. for z G X it is one 
and zero otherwise. Thus A does not depend on the momentum p. Then the scalar product in 
(|2|) reduces to 

(AV n ,y n ) = [ d^(z)\^ n (z)\ 2 , (12) 
Jx 

i. e. the probability to find the particle in X . The invariant volume element on T> is designated 
as dfi(z). With this choice of A the average a a is given by 

_ Area(X) 
° A = Avea(M) ' 

where Area(A') is the area of the domain X and Area(A^) is the area of the desymmetrized 
domain. The average 7ta is exactly what is classically expected for the mean value of the 
classical analogue of the observable A for an ergodic system. The domain X is chosen to be 
a circular domain. For the numerical integration we choose a disc with radius R having its 
center at z = 0, i. e. at the origin in the Poincare disc. The disc is then transformed by a linear 
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fractional transformation z — > (az + (3)/{(3 i 'z + a*), a,/3 € C This maps the circle centered at 
z = conformally to the domain X inside the billiards for suitably chosen a and (3. An example 
of such a domain X characterizing the operator A is shown in figure [I] for the asymmetric octagon 
as well as for the triangular billiards. These two operators are chosen for the examples shown 
in figures |2j, ||, [|, || and |TT. 



In order to understand S\(E;A) let us consider the summands in (J2J) in more detail and 
define 

d n := (A^ n ,^ n ) - a a and a n := | o„ | . (13) 

Dealing numerically with finite energies we consider the limit e —* in Sarnak's conjecture (^) 
and introduce 

c n := dnEU 4 and c n := \c n \ . (14) 

Thus the mean behaviour of a n is given by c i? n 1//4 with c being the mean value of c n . For the 
domains X shown in figure |l| we average the a n 's with a triangular shaped window function, 
denoted by (■••), with a base width of An = 100 and fitted them to f(E) = CfE^ 1 ^ with ct 
being a fit parameter. The result is shown in figure ^ for all three systems. One observes a good 
agreement in favor of Sarnak's conjecture and notes that the lower bound (||) is safely fulfilled. 

The first 3000 numbers c n , which are expected to be pseudo-random numbers with fluctu- 
ations independent of n, are shown in figure || in the case of the asymmetric octagon. One 
observes that the amplitude is of the same order over the whole range. Furthermore, no large 
exceptions are observed in agreement with Sarnak's conjecture. Also for billiard A and B no 
such exceptions are discovered. 

Numerically one observes a Gaussian distribution of the c n 's. For finite energies this distri- 
bution has, however, a non-zero mean. This is due to the imposed boundary conditions for the 
desymmetrized systems. In the desymmetrized asymmetric octagon the elliptic points, denoted 
in figure |l] as dots, play a special role. Because the hyperbolic octagons are Riemannian surfaces 
without boundary, there are no Neumann or Dirichlet boundary conditions which influence the 
eigenstates as in the case of the triangular billiards. It turns out that for the negative parity 
class, the eigenstates must vanish at the elliptic points, whereas in the case of positive parity, 
the eigenstates possess local extrema there. In the case of the regular octagon we consider as 
a special desymmetrized class the triangular billiard discussed in section ||. Here Dirichlet and 
Neumann boundary conditions influence the eigenstates significantly near the boundaries. The 
special role of the boundary influences can be very clearly seen by considering the sum over the 
modulus of the eigenstates 

W(z,E) := \*n(z)\ 2 • (15) 

E n <E 

In figure § this quantity is shown for the positive parity class of the asymmetric octagon, the 
triangular billiard A with boundary conditions incompatible with the regular octagon and for 
the triangular billiard B derived from the regular octagon. One observes deviations from an 
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equidistribution enforced by the imposed boundary conditions near the elliptic points and near 
the boundary, respectively. At the Neumann boundary the amplitude is twice the mean value 
whereas at the Dirichlet boundary the amplitude vanishes, of course. For more details, see 
The crucial point is that these boundary effects slightly modify the local expectation value of 
l^^z)! 2 far away from the boundary. This is the reason for the observed non-zero mean in 
the distribution of the c n . Alternatively, this shift can also be taken into account by using 



orbit corrections to a a in equation (g) as discussed in section |Vj. With increasing energy these 
modification decrease but since one is numerically restricted to finite energies, we take the shift 
in the Gaussian into account. For the asymmetric octagon the cumulative distribution of c n is 
shown in figure |5| in comparison with 

I f (c) = hrfc(-(c-5)/V2 l 3) (16) 

where (3 is the standard deviation and 5 the mean value of the distribution of c n . A fit excluding 
the first 50 c n 's yields the values (3 = 0.0957. . . and 5 = —0.0029. . . which is very small in 
comparison with (3. A Kolmogorov-Smirnov test gives a significance level of 55%. Neglecting 
the shift 5 leads to a significance level being almost zero. In the case of billiard A one obtains 
for the domain shown in figure |l|/3 = 0.5154 ... and <5 = -0.0277 . . . again with a significance 
level of 55%. Billiard B gives (3 = 0.4333 . . . and 5 = —0.0237 . . . with a significance level of 
94%. The Gaussian behaviour of c n implies for c n the distribution density 

1 



P(c) = \ e + e \ for c > . (17) 

Thus (3 and 5 determine the mean value c, i. e. 

2 . < . ../ 6 . /2 



e W + 6 erf -=- ~ J - for 1 5 1 < . (18) 
<t VV2/?/ V vr 



i 



Furthermore, let us define 

z n := a n -cE n 4 (19) 

which can be considered as a random variable with zero mean. The standard deviation 7 of the 
distribution of z n E\!^ = c n — c is then given by 



7 = v / /32 + 5 2_ s 2 ^ \jl--l3 for \8\<g.0 . (20) 



With the definitions a n := cE n 4 and (|19|), one has, replacing N(E) in (y) by N(E), 

s ^ ~-m (£* + • (21> 
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The first sum describes the mean behaviour S\(E; A) whereas the second sum fluctuates around 
zero. 

The mean behaviour Sk(E; A) is determined by Weyl's law 

- ^ AreaCM)^ ^ 

4-7T 

neglecting the circumference term in the case of billiards A and B, leading to 

5i(£7;A) = V a n ~ — . (23) 

Now let us consider the deviations from this mean behaviour described by the second sum 
in (^l]). Assuming that the independent random variables, the variance V(E) of the 

sum Yl,E n <E z n is given by V(E) = Y^EnKE^n with v n = jEn 1 ^ 4 . The variance V(E) can be 
estimated using Weyl's law (|22| ) 

V (E) = £ f& a W . (24) 

Thus ^(-E/jA) computed directly from the wave functions should lie within the range 

Sr dc \E;A) := ^E~\ ± =g - , (25) 
6 y N(E) Ei 

if the deviations are compared with one standard deviation and if the above assumptions are 
correct. This curve is determined from the two parameter c and 7 which in turn are determined 
by (3 and 5. A comparison of (|25"| ) with the directly computed S\(E; A) provides in this way also 
a test of the exponent | in (y)7 Figure [6| demonstrates in the case of the asymmetric octagon 
that the rate of quantum ergodicity for S\(E; A) is indeed given by the power j. In contrast to 
the discussion above, the shown model curve contains the first k true a n 's, i. e. 



Or, 

E n <E k E k <E n <E 

I v 4c 



E„<Ek 

This modification was necessary since for a given domain ^ the lowest wave functions have 
not enough structure within X in order to provide pseudo-random numbers z n . The value of k 
depends on the size of the domain X . In this case k = 50 was chosen but already k = 10 gives 
a reasonable approximation. Note that 3000 wave functions are included. 
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Alternatively, we study the rate of decline of Si(E; A) by fitting to S\(E; A) the function 



f(E) := aE-P (27) 

and compare p with the expected 7. For our three systems we choose a sequence of ten domains 
Xi. Ten circles centered at the origin of the Poincare disc are chosen which are mapped by the 
same linear fractional transformation such that the considered domains do not contain elliptic 
points. Within a sequence the radius of the domains are chosen such that the area of the 
domains increases in constant steps of 1/10-th of the area of the largest domain. The two 
sequences are shown in figure for the asymmetric octagon and for the triangular billiards. In 
the case of the asymmetric octagon the selected domains Xi have no points which are mapped 
by the transformation z — > — z to points within the same domain thus eliminating dependencies 
by the involution symmetry. In figure I the obtained Si(£;vl) are shown together with the fit 
(p7|). The functions ( p7| ) with the fitted parameters are only shown over the interval used for 
the fit. The values of the fit parameter a and p are given in table |l|. One observes for the values 
of p a good agreement with the conjectured value 4. This result clearly shows that the upper 
bound @ is valid for the energies considered here. The agreement is better in the case of the 
asymmetric octagon. In the case of the two triangular billiards slightly larger deviations occur. 
In table || the values of the fit parameter p are shown depending on the interval for which the 
fit is carried out. A tendency is observed that p approaches the conjectured value \ for higher 
energies. Using domains Xi which are centered at other locations in Ai leads to other values of 
p, some above 7 and some below. It turns out that the average over these p's is very close to | 
as it is demonstrated for ten domains having different centers in fli5| ]. Thus numerical evidence 
is given in favor of the conjecture. 

The value of a depends on the area as shown in figure ^. In order to compare the three 
systems we plotted a versus the index i = 1, . . . , 10 of the domains X{. Since Area(A^) increases 
linearly with i, a is also plotted versus the area. One observes a linear behaviour of a over a 
wide range. The asymmetric octagon shows the linear behaviour up to i = 7, and the curve 
belonging to billiard B flattens at i ~ 9. In the limit Aiea(X) — > Area(TW) the values of a n 
decrease rapidly because of the normalization of the wave functions with respect to M. This in 
turn leads to a decline of a for areas X approaching M . Thus one has three different regions for 
a, one linear region for small domains, a transition region and finally a region with decreasing 
a. 



IV The orbit theory 

In this section we consider the corrections to ~a a which are provided by the orbit theory. Fur- 
thermore, the question is posed for possible differences in the properties of the wave functions 
between arithmetic and non-arithmetic systems. This question is enforced by the fact that the 



properties of the energy level statistics are quite different [EU 22]. This fingerprint is caused by 
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asymmetric octagon 


billiard A 


billiard B 




Area 




P 






Area 


OL 


P 






Area 


OL 


P 






0.2691 


0.036 


0.251 






0.0031 


0.101 


0.234 






0.0031 


0.076 


0.211 






0.5383 


0.057 


0.251 






0.0063 


0.122 


0.209 






0.0063 


0.134 


0.221 






0.8074 


0.069 


0.249 






0.0094 


0.168 


0.215 






0.0094 


0.151 


0.212 






1.0766 


0.082 


0.250 






0.0126 


0.223 


0.224 






0.0126 


0.182 


0.215 






1.3457 


0.096 


0.257 






0.0157 


0.266 


0.226 






0.0157 


0.219 


0.220 






1.6149 


0.108 


0.262 






0.0189 


0.309 


0.229 






0.0189 


0.272 


0.229 






1.8840 


0.121 


0.268 






0.0220 


0.360 


0.234 






0.0220 


0.304 


0.232 






2.1532 


0.119 


0.261 






0.0252 


0.419 


0.239 






0.0252 


0.339 


0.235 






2.4223 


0.115 


0.254 






0.0283 


0.480 


0.244 






0.0283 


0.353 


0.234 






2.6915 


0.115 


0.253 






0.0315 


0.550 


0.249 






0.0315 


0.367 


0.233 





Table 1: The fit parameters a and p are shown for the three systems for the sequences of 
domains Xi shown in figure [?]. 

the exponentially degenerated length spectrum of the periodic orbits in the case of arithmetic 
systems. It is thus important to look at the description of the wave functions in terms of classical 
orbits. 

The orbit theory derived in pjl] allows to express a nearly arbitrarily weighted sum over 
wave functions by a corresponding sum over classical orbits. Here, in contrast to the periodic 
orbit theory regarding the eigenvalue spectrum, all non-periodic orbits are necessary, since it is 
not a trace formula. For compact Riemannian surfaces the orbit formula can be stated as 

oo ^ 

HPn) K(*) *n(z') = — £ h (coshd(z, b{z'))) , (28) 
n=o ber 

where h(p) is an even function, which is holomorphic in the strip \^sp\ < ^ +e, e > and of order 
h(p) = 0(p~ 2 ~ s ), 5 > for \p\ — ► oo. The sum on the right-hand-side runs over the Fuchsian 
group r which defines the tessellation of the Poincare disc T>. Furthermore, h(y) is the Mehler 
transform of h(p) 

poo 

h{y) = / dppt&nh(irp)P_ l+l (y)h(p) (29) 
Jo 2 

with P_i +ip (y),y > 1, being the Legendre function of the first kind. The length d(z,b(z')) is 
the length of the path connecting z and z' with the topology according to b £ T which is non- 
periodic for z = z' in general. Only in the special case in which z = z' lies on the periodic orbit 
belonging to b, the length of a periodic orbit is obtained. It is the lengths of these periodic orbits 
which are so peculiar degenerated in arithmetic systems. However, since for the description of 
the wave functions according to ( |28| ) the distribution of the lengths d(z, b(z')) is important, the 
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billiard A 


billiard B 




Area 


Pi 


P2 


P3 






Area 


Pi 


P2 


P3 






0.0031 


0.218 


0.236 


0.234 






0.0031 


0.202 


0.222 


0.225 






0.0063 


0.177 


0.206 


0.235 






0.0063 


0.234 


0.229 


0.235 






0.0094 


0.163 


0.220 


0.243 






0.0094 


0.237 


0.218 


0.221 






0.0126 


0.160 


0.229 


0.259 






0.0126 


0.233 


0.227 


0.219 






0.0157 


0.154 


0.230 


0.273 






0.0157 


0.229 


0.229 


0.227 






0.0189 


0.150 


0.237 


0.267 






0.0189 


0.230 


0.223 


0.245 






0.0220 


0.141 


0.244 


0.255 






0.0220 


0.240 


0.217 


0.251 






0.0252 


0.146 


0.252 


0.254 






0.0252 


0.250 


0.217 


0.255 






0.0283 


0.151 


0.258 


0.254 






0.0283 


0.264 


0.213 


0.249 






0.0315 


0.150 


0.259 


0.260 






0.0315 


0.276 


0.211 


0.244 





Table 2: The fit parameter p is shown for the two billiards in dependence on the fit interval: 
pi is fitted over n = 200 . . . 800, p 2 over n = 800 . . . 1400 and p 3 over n = 1400 . . . 2000. 



crucial question is, whether these lengths are, for general z and z' , differently distributed in 
arithmetic and non-arithmetic systems. 

IV.l The distribution of d(z, b(z)) 

To emphasize the role of d(z,b(z)) we specialize ( p8|) to a single wave function by choosing as 
in @ 

h(p) = coth(vrp) e-(P-P') 2 /^ 2 _ coth(vrp) e -^v'fl^ , £ >0 . (30) 

Approximating the Mehler transform for large p one gets the semiclassical expression (eq.(37) 
in H|) 

£ M P„) |*„(,)| 2 = + £ ^fff^li) e -(W^ , (31 , 

n=o v4vr fcgr , v^^smhTb 

where tj, := d(z,b(z)), and V is the Fuchsian group without the identity element which con- 
tribution has to be considered separately leading to the first term on the right side. The sum 
over r' is absolutely convergent for e > 0. Choosing e sufficiently small such that on the left- 
hand-side only the summand belonging to p' = p n contributes, allows the computation of a wave 
function solely from the Fuchsian group V. Therefore, the modulus of a wave function |5' n (z)| 
is determined solely by the lengths t^. 

Thus to address the question whether there are differences between arithmetic and non- 
arithmetic systems with respect to their wave functions, one can also investigate the properties 
of the lengths Tb of arithmetic and non- arithmetic systems. The simplest statistic is the nearest- 
neighbour length spacing distribution P(A£,z). Since the considered systems are invariant 
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under time reversal, all lengths occur twice in the sum over T which is expressed by the fact that 
Tb = r b -i with = 1. Thus in the statistic only one of the two lengths has to be included. 
The next step is the unfolding of the length spectrum using the asymptotic law 



N(z,z',t) ~ N(z,z',t) = ^e T , t^oo (32) 

where N(z,z',t) is the staircase function 

N(z,z',t) := #{b\b eT,d(z,b(z')) < r} . (33) 

The unfolded lengths are denoted by £(z, z') = N(z, z' , T b ). For the asymmetric and the regular 
octagon all elements b E T are computed which can be represented as a product of the generators 
of the group with up to 12 factors. Both groups contain only hyperbolic elements. (Elliptic and 
reflection elements occur later, when we consider the group for the desymmetrized systems.) It 
turns out that the length spacing distribution P(A£, z) is for almost all z conform to a Poisson 
distribution Pp i sson (A£) = e~^. In the case of the asymmetric octagon one finds strong 
deviations from the Poisson distribution at positions of the elliptic points of the desymmetrized 
system, which are marked in figure [l] as dots. In order to test the agreement with the Poisson 
distribution we apply the Kolmogorov-Smirnov test to the distributions in dependence of the 



position z. The results are shown in figure 1C. Here the reverse of the significance level V is 



plotted, i. e. the low regions are the domains with a Poisson distributed length spacing. It is 
clearly seen that in the case of the asymmetric octagon only the elliptic points show strong non- 
Poisson behaviour whereas the regular octagon displays also deviations along symmetry lines 
due to the reflection elements which tessellate the regular octagon in just the domains of the 
triangular billiard. However, inside the triangular billiard the lengths r& are Poisson distributed 
within the same significance level as in the asymmetric octagon excluding the elliptic points. 

Recalling that the main difference between arithmetic and non-arithmetic systems arises 
from the exponentially increasing degenerations among the periodic orbits with respect to their 
lengths, one now sees that such a behaviour does not occur with respect to r& for almost all z. 
Thus, no difference is expected between the wave functions of arithmetic and non-arithmetic 
systems. There might be, however, the possibility that the arithmetic structure of the group 
leads to higher correlations among the r^'s which are independent of z and simultaneously affect 
the sum on the right-hand-side of eq.(|3~l]), but this seems to be unlikely. 

IV. 2 Corrections to a a 

The other important application of the orbit theory is the derivation of corrections to the leading 
term a a in terms of the groups V. Let us now consider the groups for the desymmetrized 

systems. In the case of the asymmetric octagon the group is extended by S := ( ^ . j being 



a realization of the parity operation z — > —z [24]. This group contains also elliptic elements. The 
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triangular billiard is described by the reflection group which is generated by the three reflection 
elements which realize the reflections along the three edges of the billiard. This group contains 
in addition to the elliptic elements also reflection elements. (For more details, see e.g. |25|| .) 
The formula eq.(|28|) is also valid for these groups if one takes for T the above groups and takes 
the character x(b) of the group elements b E V on the right-hand-side of ( |2"8| ) into account. 
To obtain the corrections to a a we choose in eq.(^) 

hip) = e( Pmax -\ P \) , (34) 

where is the Heaviside step function. This leads to 

S~] \^ n (z)\ 2 = — dpptanh(7rp) 

+ 2^ x< ^> l dpptanh(vrp)P_i + . p (cosh(r b )) , (35) 

with N such that p^ < p max < Pn+i- The left-hand-side is now rewritten as a Riemann-Stieltjes 
integral 

N rPmax /*Pmax 

V|f„(z)f= / \V n (z)\ 2 dN(p) ~ / \V n (z)\ 2 dN(p) , (36) 

where in the last step the approximation 

dN(p) ~ dN{p) = ( ^^p + — J dp 

has been used together with Weyl's law including the circumference term. The length C is defined 
as C := C + — £~ whereby C + and C~ are the lengths of the boundary on which Neumann and 
Dirichlet boundary conditions are imposed, respectively. The comparison of the integrand of 
| with © yields 



—r { 1 + E P 

Pn I ber> 



' *"< 2) 12 " jJ$)?±C < 1 + ^ *<» "-J*^ ) 1- • (37) 

In contrast to ([H|) this expression is at most conditionally convergent. Integrating equation 
( pTD over a domain A' leads to an orbit expression for (A^! n ,^ n ) as defined in (p^). The term 
corresponding to the identity element in eq,(|37|) yields a a- The influence of the boundary, 
i. e. of the reflection elements, and of the corners, i. e. of the elliptic elements, vanishes in the 
semiclassical limit such that the properties of the a n 's are dominated by the contribution of the 
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hyperbolic elements. Thus one possible extrapolation to the semiclassical limit is provided by 
introducing 

^" = {** + A^(M) g^> L *«*-**.(«-"<'»»} (l- 2p J e&{M) ) 

(38) 

The definition of 0, C T' is arbitrarily and depends on the effects one wishes to incorporate 
in cf™ rr . In the case of the two triangular billiards A and B only the three reflection elements 
which realize the reflections at the boundary of the billiards are included in our consideration. It 
turns out that they cause the main contribution and that the contribution of further reflection 
elements as well as of elliptic elements are negligible for the domains considered here. However, 
choosing domains which contain elliptic points leads to important contributions from just the 
corresponding elliptic elements. In the case of the asymmetric octagon one has no reflection 
elements such that the elliptic elements have to be considered. The next class are the hyperbolic 
group elements whose inclusion in O is a matter of point of view. Since the orbits spread 
uniformly over the phase space in ergodic systems one suspects that they are responsible for the 
quantum ergodicity of the wave functions and thus should not be included in ( |38|) since it is just 
this effect one wants to study. 

The effect of the reflection elements is well demonstrated by the quantity 

n n := {A ^ n) (39) 
OA 

and its average («;„), where (. . .) denotes the averaging with the same triangular shaped window 
function as above. According to @ this tends towards one. The result is shown in figure |ll] 
for the asymmetric octagon and for both billiards using the domains X shown in figure [l|. As 
expected, the deviations from one are the larger the smaller n is, showing directly the influence 
of the boundary. In the case of the asymmetric octagon the deviations are negligible for n > 700 
showing only a minute influence of the elliptic elements. The larger deviations of the two 
billiards are well explained by the contribution of the reflection elements. Integrating ( |37| ) over 
the domain X and taking only the reflection elements into account provides the dashed curves 



shown in figure 11. Not only the mean behaviour of the deviations are explained but also the 



long-range oscillations. Subtracting this contribution from the K n and averaging the result gives 



the dotted curves in figure |11|. These fluctuate around one already for very small values of n. 
The remaining fluctuations are due to the contributions of other group elements not included in 
C r'. The order of the amplitude of these fluctuations has to be compared with the standard 
deviation An of the fluctuations of the K n 's which are not averaged. For billiard A and B one 
gets Ak = 0.127 and Ak = 0.106, respectively, which is orders of magnitudes larger than the 
range shown in figure 11. Computing c^ orr := \ {A^ n ,^ n ) — a^^lE 1 ^ and applying the fit (|l 



as in section III gives for the c™ rr 's for billiard A f3 = 0.5151 . . . and 5 = —0.0127. . . and for 



billiard B [3 = 0.4333 . . . and 5 = 0.0018 . . ., respectively. Thus the width of the distribution is 
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unchanged but the shift 5 of the Gaussian is drastically reduced. The significance levels turn 
out to be unaltered. 

We have computed Si(E; A) with these corrected a n 's but found no significant improvement 
concerning the fit values of p. It thus seems, that for some other reasons the two billiards 
approach slower the suggested value i than the asymmetric octagon. 



V Summary and discussion 

In this paper we have investigated the behaviour of the wave functions of three strongly chaotic 
(Anosov) systems. The rate of quantum ergodicity is studied by the limit E — > oo of S\(E\ A) 
defined in (0) and (12) and is found to be consistent with the conjectured decline proportional 
to £rV* The average of the individual summands a n , occurring in the sum in is shown 
to decline also proportional to E^ 1 / 4 . The c n := a n E\[^ , a n defined in eq.(13), are distributed 
according to a Gaussian with a non-zero mean for finite energies. Addressing the unique quantum 
ergodicity hypothesis we find no significant exceptional behaviour in individual c n 's, whereby 
we choose domains far away from elliptic points and boundary lines of the considered systems. 
These regions play a special role as it is demonstrated by the sum (fO^) of the modulus of the 
wave functions as seen in figure Sequences of domains Xi with increasing area are constructed 
for which the computed S\(E;A) are compared with f(E) = aE~ p . The values of the fit 
parameter p are consistent with the expected value \. The dependence of a on Area(Afj) is 
linear for Area(A'j) <C Area(A4) and declines thereafter towards zero for Area(A'i) — ► Area(TW) 
after passing a transition region. 

The above analysis shows no differences between the arithmetic and the non-arithmetic 
systems. The orbit theory, which allows to express the wave functions in terms of classical 
orbits, is applied to individual wave functions showing that the modulus of the wave functions 
depends only on the hyperbolic distances t& := d(z,b(z)), see eq.([3l|). Since the peculiarities 
of the statistics of the quantal levels of arithmetic systems are traced back to the exponential 
degeneration of the lengths of the periodic orbits, one has to pose the analogous question with 
respect to the wave functions. In this case it is, however, the statistical properties of the 
lengths Tb which determine possible peculiarities in arithmetic systems. The nearest-neighbour 
length spacing P(A£, z) of the unfolded r^'s is Poisson distributed within the same significance 
level in all three systems. Only at the elliptic points of the asymmetric octagon and along 
the boundaries of the billiards the distribution is non-Poisson, as the Kolmogorov-Smirnov test 



presented in figure 1C shows. Thus no peculiarities for the wave functions should be expected 
in arithmetic systems in contrast to the quantal levels. 

The deviations of (A^ n ,^f n ) from a a due to the influences of the boundary conditions are 
computed in terms of the reflection elements in the case of the two billiards. The result presented 



in figure |ll] shows that the main contribution arises from the three reflection elements realizing 



the reflection at the boundaries of the triangular billiards. 



14 



Therefore, our analysis gives strong numerical support for the quantum unique ergodicity 
hypothesis in the case of dynamical systems defined on hyperbolic surfaces independent of arith- 
metic or non-arithmetic properties. 
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Figure captions: 



Figure 1: On the left side the fundamental domain A4 of the asymmetric octagon is shown on 
the Poincare disc T>. The involution symmetry is clearly visible. The chosen operator A acts 
non-vanishingly on the grey disc X. The dots denote the position of the elliptic points on this 
surface. On the right side the triangular billiard is presented with the two boundary conditions 
chosen here. 

Figure 2: The average (a n ) over 100 neighbouring values for the domains X shown in figure |]. 
The corresponding fits f(E) = CfE" 1 ^ are shown as dashed curves. 



Figure 3: The normalized pseudo-random numbers c n are displayed for the asymmetric octagon. 



Figure 4: W(z, E) is shown for the asymmetric octagon, for billiard A and for billiard B. 
The triangular billiards are rotated such that one looks on the Dirichlet boundary. For the 
asymmetric octagon we have E = 6 000 and for the two billiards E = 200 000. 



Figure 5: The cumulative distribution I{c ) of the normalized numbers c n is shown in comparison 
with a cumulative Gaussian (dashed curve) for the same operator A as in figure |3[ 



Figure 6: Si(E;A) is shown for the asymmetric octagon as a full curve. The dashed curve is 
obtained by evaluating the sum with the mean values a n depending only on c. However, the 
first 50 correct terms a n have been included, see eq.(|26|). The dotted curves correspond to one 
standard deviation from the mean behaviour thus demonstrating the decline with the power j. 
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Figure 7: The chosen sequences of domains Xi are shown for the considered systems. 



Figure 8: S\(E; A) and the fit (27) are shown for the asymmetric octagon in a), for billiard A 
in b) and for billiard B in c). The 10 different curves belong to the 10 domains shown in figure 

0- 



Figure 9: The fit parameter a is plotted for the three systems in dependence on the area of the 
domain Xi. The full curve corresponds to the asymmetric octagon, the dotted curve to billiard 
A and the dashed curve to billiard B. 



Figure 10: The reverse of the significance level V that the unfolded length spacing of d(z, b(z)) 
is distributed according to a Poisson distribution, is shown. 



Figure 11: The quantity (/t n ) is shown for the asymmetric octagon (dashed-dotted curve) and 
for both billiards (full curves) for the domains X shown in figure (p. The lowest curve belongs to 
billiard B, and the curve belonging to billiard A lies above it. The contribution of the identity 
and the three most important reflection elements according to (|37]) is displayed as dashed curves. 
Subtracting this contribution from n n and averaging the result yields the dotted curves. 
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